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ABSTRACT 



Steam-condensation heat-transfer measurements were made 
using a 5-tube in-line test condenser with an additional 
perforated tube to simulate up to 30 active tubes. Results 
were obtained for smooth tubes, wire-wrapped tubes and 
dropwise-coated tubes. The average outside heat-transfer 
coefficient for 30 smooth tubes was 0.64 times the Nusselt 
coefficient for the first tube. A total of eight wire- 
diameter and wire-pitch combinations were tested: 1.6-mm- 

diameter wire wrapped at 16 mm, 7.6 mm and 4 mm wire pitches, 
1 . 0-mm-diameter wire wrapped at 8 mm, 4 mm and 2 mm wire 
pitches, and 0 . 5-mm-diameter wire wrapped at 4 mm and 2 mm 
wire pitches. The best bundle performance was obtained when 
the tubes were wrapped with 1 . 0-mm-diameter wire at a wire 
pitch of 4 mm. This combination resulted in an average out- 
side heat-transfer coefficient for 30 tubes that was 1.15 
times the value computed for the first tube using the Nusselt 
theory. The average outside heat-transfer coefficient for 
the 30 dropwise-coated tubes was 1.1 times the value of the 
heat-transfer coefficient for the first tube in the tube 
bundle. Utilizing either wire-wrapped tubes or dropwise- 
coated tubes, it is possible to significantly reduce the 
condenser surface area and overall size. 
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NOMENCLATURE 



2 

A q Outside heat-transfer area of one tube (m ) 

2 

A Mean Vapor flow area (m ) 

r 

2 

A^ Inside heat-transfer area of one tube (m ) 

B Coefficient defined in equation (2.8) 

0 14 

C c Correction factor (u /u ) ' 

f c w 

Cj. Calculated correction factor 

f ,c 

C. Sieder-Tate coefficient 

l 

Cp c Specific heat of water (kj/kg-K) 

Dp Average flow dimension (m) 

Inner diameter of the tube 
D q Outer diameter of the tube 

DT Temperature difference (°C) 

F Dimensionless quantity defined in equation (2.8) 

2 

g Acceleration of gravity (9.81 m/s ) 

h^ Latent heat of vaporization (kJ/kg) 

2 

h^ Inside heat-transfer coefficient (W/m K) 

h Local, outside heat-transfer coefficient for the Nth 

tube (W/m K) 

h Heat-transfer 2 coef f icient calculated from the Nusselt 

equation (W/m K) 

2 

h Q Outside heat-transf er coefficient (W/m K) 

h. Outside heat-transfer coefficient for the first 

tube (W/m 2 K) 

k c Thermal conductivity of the cooling water (W/mK) 

kp Thermal conductivity of the condensate film (W/mK) 
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Thermal conductivity of Titanium (W/mK) 

L Condensing length (m) 

LMTD Logarithmic Mean Temperature Difference (K) 
m c Mass flow rate of cooling water (kg/s) 

n Exponent defined in equation (2.8) 

N The number of tubes in a column or the tube number 

of a given tube 

Nu Water-side Nusselt number 

P^ Longitudinal pitch of the tube bundle 

P Prandtl number 

r 

P T Transverse pitch of the tube bundle 

2 

q Heat flux based on the outside area (W/m ) 

Q Heat-transfer rate (W) 

Re Water-side Reynolds number 

Re 2 ^ Two-phase Reynolds number (p^ V^D^/y^) 

R, Wall thermal resistance based on the outside area 

W (m2 K /W) 

s Inundation exponent defined in equation (5.1) 

T^ Average cooling water bulk temperature (°C) 

T ci Cooling water inlet temperature (°C) 

T cq Cooling water outlet temperature (°C) 

T c on Condensate film temperature (°C) 

Tf Average condensate film temperature (°C) 

T,. Calculated condensate film temperature (°C) 

t / c 

T Saturation temperature of steam (°C) 

S3 t 

T va p Vapor temperature (°C) 

T w Wall temperature (°C) 
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2 

U Q Outside heat-transf er coefficient (m K/W) 

V c Cooling water velocity (m/s) 

V v Vapor velocity (m/s) 

Y Dimensionless quantity defined in equation (4.21) 

GREEK SYMBOLS 

2 

y c Dynamic viscosity of the cooling water (N-s/m ) 

2 

y^ Dynamic viscosity of the condensate film (N*s/m ) 

2 

y Dynamic viscosity of water (N*s/m ) 

3 

P c Density of the cooling water (kg/m ) 

Density of the condensate film (kg/m^) 

3 

P v Density of the vapor (kg/m ) 
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I . 



INTRODUCTION 



A. HISTORICAL BACKGROUND 

Considerable interest has been generated in reducing the 
size and the weight of propulsion systems for naval applica- 
tions. Advances in condenser design could do much to reduce 
the size and the weight of the propulsion plant. Measures 
to raise the condensing-side heat-transfer coefficient (of 
condenser tubes) is one way to achieve this reduction in 
condenser size. This reduction, however, comes at a price. 
This is usually due to an increase in the pumping power or 
due to an increase in the initial cost of the tubes. For 
naval applications, where the size of a vessel may depend 
upon the size of the condenser (in a submarine, for example) 
this reduction in the size is justified, even at the greater 
cost . 

Search [Ref. 1], at the Naval Postgraduate School, 
investigated the present condenser design process, and 
examined the potential benefits that might occur if heat- 
transfer enhancement was established in the condenser. He 
concluded that reductions of as much as forty percent in the 
size and weight of condensers are possible. This is depen- 
dent, of course, on the heat-transfer-enhancement technique 
utilized. Much further research work at the Naval Post- 
graduate School has been directed toward these heat-transfer 
enhancement techniques. 
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Beck [Ref. 2], Pence [Ref. 3], Reilly [Ref. 4], Fenner 
[Ref. 5], and Ciftci [Ref. 6] conducted research employing 
a single- tube test condenser. Their research concluded that 
the overall heat- transfer coefficient of enhanced tubes may 
be as much as twice that for smooth tubes of similr geometry. 
In a separate report, Marto, Reilly and Fenner [Ref. 7] re- 
ported that most of the increase in the overall heat-transfer 
coefficient was on the cooling-water side and was due to an 
increase of the turbulence and the swirl, as well as to an 
increase in the inside surface area. Only a small increase 
occurred on the steam side. 

Present-day steam condensers utilizing smooth tubes are 
limited in their thermal efficiency, due primarily to the 
large thermal resistances occurring on the tube side of the 
condenser. It is possible, however, by utilizing enhanced 
tubes, to increase the inside heat- transfer coefficient by 
100 percent or more. The corresponding increase in the out- 
side heat-transfer coefficient, however, is less than 50 
percent. In studying ways to further increase the outside 
heat-transfer coefficient, Webb [Ref. 8] reported that 
conduction across the condensate film is the primary thermal 
resistance in film condensation. This thermal resistance is 
usually larger than the thermal resistance of the tube wall, 
that attributed to fouling or that due to noncondensable 
gases. It is possible to reduce this thermal resistance by 
utilizing a geometry that reduces the film thickness. This 
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reduction of the thermal resistance would mean a correspond- 
ing increase in the outside heat-transfer coefficient. 

For large tube bundles, condensate inundation is present 
and must be considered when attempting to increase the 
overall heat-transfer coefficient. Thomas [Ref. 9] wrapped 
wire around smooth tubes in a helical manner and tested the 
condensation of ammonia in a large tube bundle. Increases 
of as much as 200 percent in the oustide heat-transfer 
coefficient over that predicted by Nusselt [Ref. 10] were 
measured. This increase was attributed to the effect of 
surface tension drawing the condensate to the wire and 
acting as a condensate run-off channel. In a paper by 
Cunningham [Ref. 11], "roped" tubes were considered and 
again the increase in the outside heat-transf er coefficient 
was attributed to condensate drainage, while an increase in 
the inside heat-transfer coefficient was attributed to the 
increase in the inside convective coefficient, due to the 
increased turbulence. Kanakis [Ref. 12] tested both smooth 
and roped tubes, with and without a wire wrap in an in-line 
tube bundle simulating up to 30 tubes. Adding the wire wrap 
on the smooth tube increased the average, outside heat- 
transfer coefficient for the 30 tube bundle by 50 percent, 
while adding the wire warp to the roped tube increased the 
average outside heat-transfer coefficient for the 30 tube 
bundle by more than 35 percent. 

For condenser tubes, the increase in the outside heat- 
transfer coefficient may also be accomplished by promoting 
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dropwise condensation. Tanasawa {Ref. 13] reviewed dropwise 
condensation and discussed the methods for promoting these 
dropwise conditions. He concluded that the use of organic 
polymers (such as Teflon) was the most promising of all the 
dropwise-promoting techniques. Investigations of organic 
coatings by Brown [Ref. 14] have shown that enhancements of 
up to 180 percent are possible. Of primary concern, however, 
is that these coatings have very low conductivities and must 
therefore be ultra-thin. Another concern is that these 
coatings must be strongly adherent and sufficiently tough to 
withstand the conditions of their assembly and their use. 
Holden [Ref. 15] investigated numerous dropwise-promoting 
coatings. His tests were based not only on the ability of 
the coating to promote dropwise condensation, but also on 
the ability of the coating to sustain this dropwise perfor- 
mance over an extended period of time. In addition, the 
coatings that were able to sustain this dropwise performance 
were evaluated for their heat-transfer performance. His 
results indicated that the outside heat-transfer coefficient 
for a single tube could be increased by a factor of from 
five to eight through the use of polymer coatings. 

Since the use of a wire wrap can increase the heat- 
transfer performance of a tube bundle, an important question 
arises: Is there some optimal wire diameter and wire pitch 

combination? In addition, during dropwise conditions what 
is the effect of condensate inundation upon heat transfer in 
a tube bundle? 
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B. OBJECTIVES 



The objectives of this thesis were therefore to: 

1. Conduct steam-condensation tests to determine the 
steam-side heat-transf er coefficient for 16-mm-o.d. 
smooth titanium tubes, 

2. Confirm the heat-transfer-performance measurements 
of Kanakis [Ref. 12] on 16-mm-o.d. smooth titanium 
tubes with a 1.6-mm-o.d. titanium wire wrapped at a 
pitch of 7.6 mm, 

3. Conduct heat-transfer-performance measurements to 
determine the effect of wire pitch and wire 
diameter on the steam-side heat-transf er 
coefficient with inundation, and 

4. Conduct heat-transfer-performance measurements to 
determine the effect of inundation on a dropwise 
coated tube and compare this performance with the 
performance of a smooth tube. 
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II. THEORETICAL AND EMPIRICAL BACKGROUND 



The basis for the analysis of film condensation on a 
horizontal tube was set forth by Nusselt in 1916. His analy- 
sis was, however, for laminar film condensation on a single 
horizontal tube. Nusselt 's analysis yielded the well-known 
relationship for the heat-transfer coefficient: 



Nu 



= 0.725 



-.3 2 . 

k f P f h f q g 
D o p f (T sat- T „> 



1/4 



( 2 . 1 ) 



This relationship is subject to the following restric- 
tions, as stated by Nobbs [Ref. 16]: 

1. The wall temperature is constant, 

2. The flow is laminar in the condensate film, 

3. The film thickness is small compared to normal 
tube diameters, 

4. All fluid properties are constant within the 
condensate film, 

5. Heat transfer in the film is by conduction, 

6. All forces due to hydrostatic pressure, surface 
tension, inertia, and vapor/liquid interfacial shear 
are negligible when compared to viscous and 
gravitational forces, and 

7. The vapor/liquid interface and the surrounding steam 
are at the saturation temperature. 

Jakob [Ref. 17] extended the Nusselt analysis to film 
condensation on a vertical in-line column of horizontal 
tubes by assuming that all the condensate from a tube drains 
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as a laminar sheet onto the tube below it. This is depicted 
in Figure 2.1a. In this idealized situation, the average 
coefficient for a vertical column of N tubes was predicted 
to be : 



h = 0.725 

N 



■4 4 v 9 



1/4 



D o N >*f 



(T — T ) 

sat w 



( 2 . 2 ) 



Combining equations (2.1) and (2.2) yields the Nusselt 
idealized theory for the average coefficient compared to the 
coefficient of the top tube: 



h N /h l 



= N 



-1/4 



(2.3) 



In terms of the local heat-transfer coefficient for the Nth 
tube, this result becomes: 

h N /h 1 = N 3/4 - (N - 1 ) 3/4 (2.4) 

Realizing that condensate does not flow in a laminar sheet 
but by discrete droplets, Kern [Ref. 18] proposed a less- 
conservative relationship to account for the ripples and 
turbulence introduced into the condensate film. This is 
depicted in Figure 2.1b. The Kern relationship is: 

hl/h. = N _1/6 (2.5) 

N 1 
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a. Nusselt's Idealized 
Laminar Flow Model 



b. Kern's More Realistic 
Flow Model 




c. Eissenberg's Side-Drainage Flow Model 
Figure 2.1 Representations of Condensate Flow 
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or in terms of the local heat-transfer coefficient for the 



Nth tube: 




N 5/6 - (N - 1) 5/6 



( 2 . 6 ) 



Eissenberg [Ref. 19] did extensive experimentation to 
investigate the effects of steam velocity, condensate inun- 
dation and noncondensable gases on the condensation heat- 
transfer coefficient. He theorized that condensate does not 
always drain only onto tubes aligned vertically, but can be 
diverted sideways, especially in staggered tube bundles. 

The condensate draining sideways strikes the lower tubes on 
their sides rather than on their tops. This is depicted in 
Figure 2.1c. Since more heat is transferred from the top of 
a tube than from its bottom, the net effect of inundation is 
less severe. He predicted the following relationship: 



Extensive experimental research into the effect of 
condensate inundation has been conducted. However, the data 
exhibit a substantial amount of scatter as shown in Figure 
2.2 by the cross-hatched area. Berman [Ref. 20] conducted a 
compilation of film condensation data on bundles of horizontal 
tubes and identified the following variables as important 
in causing the scatter: 




0.60 + 0.42 N 1/4 



(2.7) 
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Figure 2.2 Theories and Data for Condensate Inundation 



1 . 



Bundle geometry (in-line or staggered) , 

2. Tube spacing, 

3. Type of condensing fluid, 

4. Operating pressure, 

5. Heat flux, and 

6. Local vapor velocity. 

Marto and Wanniarachchi [Ref. 21] noted that other factors 
can have an effect on the data such as noncondensable gases, 
the direction of the vapor flow, partial dropwise condi- 
tions, an insufficient abount of steam reaching the lower 
tubes in the bundle and difficulty in the measurement of 
the local condensing coefficient. 

Very small amounts of noncondensable gases can result in 
significant reductions in the condensation heat-transfer 
rate. Summaries of the phenomenon have been made by Chisholm 
[Ref. 22], and Webb and Wanniarchchi [Ref. 23], The consen- 
sus among them is that noncondensable gases impose an added 
thermal resistance, since the vapor molecules must diffuse 
through a gas layer prior to reaching the condensing surface. 
These noncondensable gases can also lead to regions where 
the tubes are inoperative in a condensing role. Nonconden- 
sable gases have an adverse effect on the condenser perfor- 
mance, and must be taken into account when designing a tube 
bundle . 

The motion of vapor within a condenser affects the film 
condensation process because of its effect on the surface 



23 



shear between the vapor and the film, and the resulting effect 
on vapor separation. Although the results may be different 
depending upon the orientation of the steam flow, the general 
effect of an increase in vapor velocity is a corresponding 
increase in the condensing heat-transf er coefficient. There 
exist a number of both theoretical and empirical equations 
representing vapor-shear effects on the condensate heat- 
transfer coefficient [Refs. 24, 25 and 26]. 

It is clear that the measurements that were made during 
this thesis had both inundation and vapor-shear effects. 

Marto [Ref. 27] stated that these two effects are difficult 
to separate from one another. During this study, however, 
an attempt was made to separate these two effects. For this 
purpose, data were taken on the top tube of the bundle with 
vapor velocity as a variable, and the data were correlated 
using an expression similar to a correlation suggested by 
Fuj ii [Ref . 28] . 



1/2 

Nu/Re 24> 



BF 



n 



( 2 . 8 ) 



This correlation is based on numerous data for both 
staggered and in-line tube bundles. To represent the test- 
condenser tube bundle, the constants B and n in equation 
(2.8) were computed using a least-squares technique. 
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III. EXPERIMENTAL APPARATUS 



The test facility was designed and built originally by 
Morrison [Ref. 29], and modified by Noftz [Ref. 30] and 
Kanakis [Ref. 12] to simulate an active tube column of up to 
30 in-line tubes. A detailed description of all the com- 
ponents is contained in Ref. 30 and further modifications 
are described in Ref. 12. The descriptions in this report 
are brief, with particular focus given to the modifications 
undertaken by the author. 

A. STEAM SUPPLY 

The steam supply system is shown in Figure 3.1. House 
steam flows through a supply valve (MS-3) , into a steam 
separator, and is throttled down by another valve (MS-4) 
to operating conditions. The steam then passes through an 
orifice and into the test condenser diffuser before entering 
the test condenser where it flows through a simulated in-line 
tube bundle . 

B. TEST CONDENSER 

The dimensions of the test condenser shown in Figure 3.2 
were unchanged from Noftz 's original design. The condenser 
consisted of five active tubes, twelve dummy tubes and one 
perforated tube. The tubes were positioned in the test con- 
denser with a pitch-to-diameter ratio of 1.5. The active 
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PRESSURE U-TUBE PRESSURE 

GAUGE MANOMETER GAUGE 




Figure 3.1 Schematic Diagram of Test Apparatus 
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Figure 3.2 Sketch of Test Condenser 
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condensing length of each tube was 305 mm. A double-walled 
glass viewing window was provided on the front of the test 
condenser to allow for visual observation of the condensation 
process. Although the viewing window was designed to have 
heated air flowing between the glass walls, this feature was 
not utilized since no fogging of the glass was observed. A 
modification of the test condenser was made to minimize the 
steam flow along the test condenser side walls. The modifica- 
tion consisted of the addition of twelve baffle plates, placed 
as illustrated in Figure 3.2. 

C. TEST-CONDENSER TUBES 

Ten different types of active tubes were tested in this 
set of experiments. The tubes were manufactured by Wolverine 
Division of Universal Oil Products, and were made of titanium. 
All were of 16 mm o.d. and had a minimum wall thickness of 

0.89 mm. The ten types of tubes are as listed below: 

1 . Smooth tubes , 

2. Smooth tubes wrapped with 1.6-mm-o.d. titanium wire 
at three different pitches (16 mm, 7.6 mm and 4 mm) , 

3. Smooth tubes wrapped with 1.0-mm-o.d. titanium wire 
at three different pitches (8 mm, 6 mm and 4 mm) , 

4. Smooth tubes wrapped with 0.5-mm-o.d. titanium wire 
at two different pitches (4 mm and 2 mm) , and 

5. Smooth tubes with a dropwise coating. 

The actual wrapping of the wire around the smooth tubes 
was performed by the author of this thesis. A detailed 
description of the procedure is included in Chapter III. 
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The dropwise-coated tubes were sent to General Magna- 
plate Corporation where a dropwise coating, commercially 
referred to as "Nedox," was applied. First an electro- 
deposited nickel-cobalt substrate was applied to the tube 
surface, and then this substrate was infused with a micro- 
thin Teflon layer. A heat-treating cycle was employed to 
ensure thorough infusion of the substrate with the Teflon. 

This coating is highly non-wetting, has a low coefficient 
of friction and is heat-resistant. The specifics of the 
"Nedox" process are proprietory, but the manufacturer claims 
control of the surface thickness to 0.0025 mm. For this thesis 
the requested surface thickness was 0.0075 mm. 

D. PERFORATED TUBE 

Each set of five active tubes had its own corresponding 
perforated tube. This perforated tube was located above the 
uppermost active tube; see Figure 3.2. Figure 3.3 shows a 
typical pattern for the perforations of a tube, and is an 
illustration of the supply-water tubing at the test condenser. 
This particular arrangement of the tubing was used in an 
attempt to provide an even flow from the entire length of the 
perforated tube. A heated tube was used to provide the 
water supply for the perforated tube. A separate rotameter 
was provided in order to control the rate of flow of condensate 
through the perforated tube. 
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Perforated- tube water supply flow shown by arrows 



E. CONDENSATE SYSTEM 

The condensate system remained unchanged from Kanakis' 
modifications. Figure 3.1 shows the system as designed by 
Noftz and as modified by Kanakis. As steam condenses in the 
test condenser, it runs into and collects in the test con- 
denser hotwell. The hotwell is a calibrated cylinder equipped 
with a sight glass, providing a visual measure of the 
condensate being produced. Measurement of the condensate 
collection rate provided the flow rate to be supplied to 
the perforated tube. 

F. COOLING WATER SYSTEM 

Although the cooling water system remained unchanged 
from Noftz original design, the five rotameters were recali- 
brated after Kanakis' data were taken. Each active tube had 
its own rotameter and regulating valve, and it was possible 
to control the coolant velocity through the active tubes up 
to values of 5 m/s. All data for this thesis, however, were 
taken at a constant coolant velocity of 1.56 m/s. The inlet 
temperature of the cooling water was maintained at a nearly 
constant temperature by the use of a large supply tank and 
a separate cooling tower to dissipate excess heat to the 
atmosphere. In order to provide a bulk temperature at the 
outlet of each active tube, a mixing chamber was installed 
just upstream of the temperature monitoring point. 
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G. INSTRUMENTATION 

1 . Flow Rates 

Rotameters were used to measure the flow rate of the 
cooling water to each active tube and to the perforated tube. 

2 . Temperatures 

Stainless-steel-sheathed copper-constantan thermo- 
couples were used as the primary temperature monitoring 
devices. Table 1 lists the locations monitored by each of 
the thermocouples. Calibration, as performed by Kanakis, 
was utilized for all data taken. 

The perforated-tube water supply was not maintainable 
at a constant temperature, nor was it quickly adjustable to 
a desired temperature. The system, as designed, contained 
a Gulton Industries, West 20 temperature controller. Although 
the manufacturer's accuracy was listed as 1.25 degrees F, 
the system response was so poor that the temperature of the 
perforated-tube supply water varied by as much as 10 degrees 
C from the beginning of the data run to the completion of 
the data run. 

3 . Pressure 

A Bourdon-tube pressure gauge was used to measure the 
house steam supply pressure, while a compound pressure gauge 
was used to measure the steam pressure downstream of the 
orifice. A pressure transducer was used to measure the 
pressure in the test condenser, and a U-tube mercury manom- 
eter was fitted to measure the pressure drop across the 
orifice . 
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TABLE 1 



Thermocouple Monitoring Locations 



LOCATION 


CHANNEL 


T . 

Cl 


#1 


000 


T - 

Cl 


#3 


001 


T . 

Cl 


#5 


002 


T 

CO 


n 


003 


T 

CO 


#1 


004 


T „ 

CO 


#1 


005 


T 

CO 


#2 


006 


T 

CO 


#2 


007 


T ^ 
CO 


#3 


008 


T 

CO 


#3 


009 


T 

CO 


#4 


010 


T 

CO 


#4 


Oil 


T 

CO 


#4 


012 


T 

CO 


#5 


013 


T „ 

CO 


ft 5 


014 


T 

sat 




015 


T 

sat 




016 


T 

con 


017 


T 

vap 


018 
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4 . Data Collection and Display 

A Hewlett-Packard model 3054A Automatic, Data- 
Acquisition System, with a HP 2671G printer was used to 
record and display the thermocouple and transducer readings. 
The pressure transducer was assigned to channel 19 of the 
data-acquisition system, with the thermocouples assigned to 
the channels as indicated in Table 1. 
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IV. EXPERIMENTAL PROCEDURES 



A. CALIBRATION 

1 . Rotameters 

Calibration was performed on the five rotameters 
that supplied the cooling water to the active tubes. This 
calibration was performed using a weighing tank and a stop- 
watch, and a least-squares-f it calibration line was generated 
for each of the rotameters. 

2 . Steam Orifice 

To calibrate the steam orifice the following steps 
were taken: 

1. Set a steam flow rate, 

2. Adjust the cooling water flow rate to condense all 
the incoming steam, 

3. Measure the pressure drop across the orifice, 

4. Measure the condensate collection rate, 

5. Repeat steps 1 through 4 for different steam flow 
rates, and 

6. Develop a least-squares-f it straight line for the 
steam mass flow rate vs. pressure drop. 

B. PREPARATION OF CONDENSER TUBES 

Prior to the winding of the titanium wire, the outside 
surface of the tubes was buffed using steel wool unitl all 
evidence of surface oxidation was removed. If any oxidation 
was evident on the inner surface of the tube, a test-tube 
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brush was used to apply a 50-percent sulfuric-acid solution 
to the inner surface of the tube. The tube was then 
thoroughly rinsed using tap water. 

The procedure followed in the wrapping of the wire onto 
the smooth tubes was the same for every wire diameter and 
pitch combination tested. The wire was first welded to the 
tube and then, under a constant tension of 22 Newtons, the 
wire was guided (using a metal plate to control the spacing 
between the wire wraps) to the proper pitch as the tube was 
manually turned. Upon completion of the wrapping, the free 
end of the wire was then welded to the tube to yield a 
helically-wrapped surface 305 mm in length. After the wire 
wrapping of the tubes was completed and just prior to the 
installation of the tubes in the test condenser, the tubes 
were cleaned inside and out using a brush and a biodegradable 
detergent. A thorough rinsing was agin undertaken, and the 
tubes were checked to ensure that they displayed proper wetting 
characteristics. If any non-wetting areas were evident, a 
50-percent solution of sodium hydroxide mixed with an equal 
amount of ethyl alcohol and heated to about 80 degrees C, was 
brushed onto the tube surface. After rinsing with tap water, 
the tube was ready for installation in the test condenser. 

Upon completion of the installation of the tubes in the test 
condenser, steam was introduced into the test condenser. 

Again the tubes were visually examined. If there were any 
areas that displayed dropwise condensation, the tubes were 
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washed from above using the perforated- tube water supply. 

This was sufficient in all cases to restore filmwise 
condensation . 

C. SYSTEM OPERATION 

The system was operated in accordance with the operating 
instructions outlined in Appendix A of Ref. 17. The system 
was deemed to be operating at a steady-state condition when 
the cooling water inlet temperature was steady. The cooling 
water inlet temperature was monitored using the HP-3054A 
Automatic Data Acquisition System, and was considered steady 
when two consecutive readings (at five-minute intervals) 
varied by less than 0.1 degrees C. 

One hour was usually required from system start-up to 
the steady-state condition. Upon reaching the steady-state 
condition, a data run was begun. Any changes in the cooling 
water flow or changes in the perforated-tube water flow were 
accompanied by a five-minute wait for system restabilization, 
prior to taking any further data. The duration of each data 
set was approximately one minute , and a completed data run 
consisted of 30 separate data sets. During each data set, 
the condensation collection rate was measured. For each 
inundation condition, five consecutive data sets were taken 
and the average values were computed. After each set of 
five consecutive data sets, the average condensate collection 
rate was used to determine the rate of flow of the water into 
the perforated tube which was to be used for the next set of 
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data (the average amount of condensate collected during a 
data set was introduced into the perforated tube for the 
next inundation condition). For example, to simulate tubes 
6 through 10, the perforated-tube water flow rate was set 
equal to the condensate collection rate for tubes 1 through 
5. A data run was considered completed upon the simulation 
of inundation conditions through the 30th tube. 

Thermocouple readings and pressure-transducer readings 
were taken automatically by the data-acquisition system, 
while the settings of the rotameters and the test condenser 
initial and final hotwell levels were entered into the 
computer using the keyboard. Initial and final hotwell 
levels were taken for a time interval of one minute. 

The following conditions were the operating conditions 
for all data taken for this thesis: 

1. The coolant velocity was 1.55 ± 0.05 m/s, 

2. The coolant inlet temperature was 24 degrees C. (±3 
degrees C. depending on the day of the run, but 
constant during a given run) , and 

3. The saturation temperature was 100.5 ± 0.5 degrees C. 
D. DATA- REDUCTION PROGRAM 

A computer program was utilized to process and plot all 
raw data. The program was written in BASIC language and 
was run on an HP-9826 computer. A listing of this program 
is included as Appendix C. 
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E. HEAT-TRANSFER-COEFFICIENT CALCULATION 



1. Inside Heat-Transfer Coefficient 



The inside heat-transfer coefficient was computed 
using the Sieder-Tate equation described in Holman [Ref. 



31] : 




C. Re 



1 




(4.1) 



where the coefficient C^ was computed using a modified 



Wilson plot [Ref. 12]. Using this technique, Kanakis [Ref. 

12] obtained a Sieder-Tate coefficient of 0.029 ± 0.001. 

The method Kanakis used, however, did not include the variation 
of the condensate film properties as a function of heat flux. 

A slightly different, modified Wilson plot, as described by 
Nobbs [Ref. 16], that does account for this variation in 
the film properties, was used in this thesis. Appendix D 
gives a description of this method. Reprocessing Kanakis' 
data with this method, the Sieder-Tate coefficient was deter- 
mined to be 0.028 ± 0.001. 

2. Outside Heat-Transfer Coefficient 



computed using an energy balance, and can be expressed in 
terms of an overall heat-transfer coefficient by: 



The heat-transfer rate to the cooling water can be 



Q 




U A LMTD 
o o 



(4.2) 



where 
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LMTD = 



T -T . 
co ci 



ln((T sat- T ci )/(T sat- T co> ) 



(4 . 2a) 



Solving for the overall heat-transfer coefficient gives: 



U 



= m 



C /A 
pc o 



In ( (T 



sat 



- T . ) / ( T 
ci' 7 v sat 



-T ) ) 
co 



(4.3) 



The overall thermal resistance can be written as the sum of 
the internal, wall and external resistances: 



w 



U A 
o o 



h . A. 
l l 



h A 
o o 



(4.4) 



where 



R = D ln(D /D . )/2 k 
w o o' i m 



(4.4a) 



or from this result, the outside heat-transfer coefficient 
may be written: 



h 



o 



^o 



D /D. 
o 1 




(4.4b) 



F. DAT A- RE DUCT ION TECHNIQUE FOR THE OUTSIDE HEAT-TRANSFER 
COEFFICIENT 

A listing of the computer program used in the calculation 
of the outside heat-transfer coefficient can be found in 
Appendix C. The steps required in the technique are outlined 
below : 
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1. Calculate the average bulk cooling water temperature: 



T. = (T + T . ) / 2 
b co ci 



(4.5) 



2. Calculate the cooling water velocity: 



V = m /A. p 
c c' 1 M c 



(4.6) 



3. Calculate the cooling water Reynolds number 



Re 



V 



D./y 

1 c 



(4.7) 



4. Calculate the heat transferred to the cooling 
water using the left-hand side of equation (4.2). 

5. Calculate the heat flux: 



q = Q/(tt D q L) 



(4.8) 



6. Since the film temperature was not a known quantity, 
an iterative scheme was employed to calcualte the 
film temperature, as indicated below: 

a) Assume a film temperature (say = T , 

b) ' Calculate the Nusselt coefficient. 



Nu 



= 0.651 



F k f P f h fq ^ 
“f D o 5 



1/3 



(4.9) 



c) Evaluate T , and 
r / c 



f ,C 



= T 



sat 2 h 



(4.10) 
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d) Repeat steps a) through c) until T f is 
approximately equal to T^. ,c 

7. Calculate the inside heat-transfer coefficient: 



h. = k C. Re 0 * 8 Pr 1//3 c,/D. ( 

1 c i f 1 

where is the correction factor (y /u ) 8 ' 3 ^. ( 

f c p w 

8. Since the correction factor is evaluated at the 
inner wall temperature, an iterative scheme was 
employed to calculate the wall temperature as 
indicated below: 

a) Assume a value for (e.g., 1.0), 

b) Calculate h^ using equation (4.11), 

c) Calculate a water-side temperature drop. 



DT 



q D./h. 

^ l l 



D 

o 



( 



d) Calculate a new correction factor, and 



'f ,c 






0.14 



e) Repeat steps a) through d) until C^ c is 
approximately equal to C^ . ' 

9. Calculate the log-mean-temperature difference using 
equation (4 . 2a) . 

10. Calculate the overall heat-transfer coefficient: 



U 



o 



4.11) 

4.11a) 



4.12) 

4.13) 



q/LMTD 



(4.14) 



11 . 



Calculate the outside heat-transfer coefficient using 
equation (4 . 4b) . 



12. Calculate the normalized, local, outside heat- 
transfer coefficient. 




13. Calculate the normalized, average, outside heat- 
transfer coefficient. 




NOTE: In the above formulation, the thermal resistances due 

to noncondensable gases and any fouling were assumed to be 
negligible . 

G. DATA- REDUCTION TECHNIQUE FOR THE VAPOR-SHEAR CORRELATION 
For the calculation of the vapor-shear correlation, 
measurements on only the first active tube in the tube 
bundle were taken. 

1. Calculate the inside heat-transfer coefficient as 
outlined above. 

2. Calculate the outside heat-transfer coefficient 
as outlined above. 

3. Calculate a temperature correction: 



DT 



q/h 



(4.15) 



o 



4. Assume a film temperature: 



T 



f 



Tgat " DT / 2 



(4.16) 
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5. Calculate the steam velocity: 



V 

v 



m v /A r 
v v f 



where 



A F = 2 L < P t P L D o /4)/P L 



Ap is the Mean Vapor Flow Area defined by Nobbs 
[Ref. 16] and depicted in Figure 4.1. 

6. Calculate the two-phase Reynolds number: 



Re 



24> 



p, V D /y_ 
f v o f 



7. Calculate the dimensionless quantity F: 



F = D h_ g/V 2 k c DT 

o H f fg v f 



8. Calculate the Nusselt number: 



Nu = h D/k 
o o / f 

9. Calculate Y: 

Y = Nu/Re 2( ^ 1/2 



(4.17) 



(4.17a) 



(4.18) 



(4.19) 



(4.20) 



(4.21) 



44 




To calculate the mean flow area the following steps 
are required: 

1) Calculate the cross-hatched area above, 

AREA = P T P T - IT D 2 / 4 
1 L o 

2) Divide by longitudinal pitch to obtain an 
average flow dimension 

Dp = AREA / P L 

3) Multiply by the condenser length and the 
number of flow paths (2 flow paths in the 
case of this thesis). 

A F - 2 L (P T P L - T D q 2 / 4) / P L 

Figure 4.1 Mean Vapor Flow Area 
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10 . 



Repeat steps 1 through 9 at different steam mass 
flow rates into the test condenser. For this thesis 
10 different steam mass flow rates were considered, 
and three measurements taken at each mass flow rate. 

11. Starting with an assumption that the data would 
be of a form similar to that postulated by Fujii 
[Ref. 28] , an expression was assumed: 



Y = B F n ; F as defined in equation (4.19) (4.22) 



12. Perform a least-squares-f it of the data and 
determine the constant B and the exponent n. 
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V. RESULTS AND DISCUSSION 



A. RESULTS BEFORE AND AFTER TEST CONDENSER MODIFICATION 

Initially, data were taken with the test condenser used 
by Kanakis [Ref. 12], with no modifications. As can be seen 
from Figure 5.1, the nondimensionalized heat-transfer coeffi- 
cient exhibits a saw-toothed variation over the inundation 
range considered. The particular wire pitch and wire diameter 
combination chosen for Figure 5.1 was representative of the 
pattern apparent in the data sets taken prior to the test- 
condenser modification. 

It was determined that the original design of the test 
condenser allowed too much steam to bypass the five active 
tubes, thereby allowing the steam to flow between the walls 
of the test condenser and the two rows of dummy tubes. Since 
steam was bypassing the active tubes, inadequate steam was 
flowing in the vicinity of the active tubes. This led to 
artificially low heat-transfer readings in the lower tubes 
of the tube bundle. To minimize this problem, baffle plates 
were installed on each side of the test condenser as des- 
cribed previously (Figure 3.2). 

After the placement of the baffle plates, the previously 
tested tube sets were retested. As can be seen in Figure 
5.1, the saw-toothed variation of the curve is much less 
pronounced and appears to be nearly corrected. The net 
effect of the removal of the saw-toothed variation by the 
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Figure 5.1 Data With and Without Baffles 



installation of the baffle plates is an increase in the 
measured heat-transfer coefficient. This, in turn, means 
that the effect of inundation is not actually as great as 
had been earlier suspected. All data represented in this 
report will, hereafter, be with the baffles in place. 



B. VAPOR-SHEAR CORRELATION 

As discussed previously, the vapor shear correlation was 
assumed to be of the form postulated by Fujii [Ref. 28] — 
see equation (4.8). Figure 5.2 is the graph of Y vs. F for 
the first active smooth titanium tube. For this tube, the 
value of B was determined to be 0.834 and the value of n was 
determined to be 0.166. These results fall about 15 percent 
below the Fujii correlation [Ref. 28] (B = 0.96 and n = 0.20). 
There are several possible explanations for this difference. 

1. Both the value of Y and the value of F are partly 
based on the steam velocity, with a higher value of 
steam velocity resulting in a lower value for both Y 
and F. During this thesis, the steam velocity was 
only varied from 1.4 to 2.1 m/s. The corresponding 
range of F was therefore from 3.6 to 1.5, while the 
corresponding range of Y was from 1.1 to 0.8. Since 
the range of both F and Y were considerably more exten- 
sive for the calculation of the Fujii correlation, 
this is a possible explanation for the discrepancy. 

2. For the Fujii correlation, both in-line and staggered 
tube bundles were considered. This could also be a 
possible reason for the discrepancy. 

3. The direction of the steam flow considered for the 
determination of the Fujii correlation included 
vertically downward flow, vertically upward flow and 
also horizontal flow. This difference in flow 
directions is another possible explanation for the 
discrepancy between the correlation determined in 
this report and the Fujii correlation. 
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Figure 5.2 Vapor-Shear Correlation 



4. As mentioned previously, baffle plates were installed 
in order to reduce the steam bypassing the active 
tubes. Even though the baffles reduced the amount of 
steam bypassing the active tubes, this is no guarantee 
that the steam still able to bypass the active tubes 
could not have had some influence on the vapor-shear 
correlation calculated in this thesis. 

Since Reference 28 does not specifically say how the 
flow area used in his correlation was determined, this is 
another possible area for discrepancy. In addition, Fujii 
[Ref. 28] does not show which tube in the bundle was used 
for the formulation of his correlation. The actual location 
would have an important effect on the vapor velocity and 
could be another possible area for discrepancies to develop. 

Vapor-shear correlations were determined for all of the 
tube sets tested for this report. Table 2 is a summary of 
the constants determined for the vapor-shear correlations. 

C. THE EFFECT OF INUNDATION ON A SMOOTH TUBE BUNDLE 

1 . With Vapor-Shear Contribution Present 

Figure 5.3 shows the variation of the normalized, 
local heat-transfer coefficient for up to 30 tubes. The 
value for the normalized, local heat-transfer coefficient 
for the 30th tube is approximately 45% above the value pre- 
dicted by Nusselt theory, but very near the value predicted 
by the Kern relationship. As discussed in Chapter II, the 
Nusselt theory for a tube bundle is based on idealized laminar 
flow from tubes above to tubes below. The idealizations of 
Nusselt result in a conservative estimate of the heat-transfer 
coefficient, so the fact that the data are 45% above the 
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TABLE 2 



TUBE 

SET 

1 

2a 

b 

c 

3a 

b 

c 

4a 

b 

5 



Vapor-Shear Correlations 



WIRE 



DIAMETER (mm) 


PITCH (mm) 


B_ 


n 


Smooth tube 




0.834 


0.185 


1.58 


16 


0.835 


0.178 


1.58 


7.6 


0.837 


0.189 


1.58 


4 


0.725 


0.192 


1.01 


8 


0.916 


0.161 


1.01 


6 


0.871 


0.173 


1.01 


4 


0.876 


0.220 


0.50 


4 


1.036 


0.200 


0.50 


2 


0.915 


0.192 


Dropwise-coated tube 


0.796 


0.166 
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Figure 5.3 Inundation Effect on a Smooth Tube Bundle 

(Normalized, Local Heat-Transfer Coefficient) 



Nusselt prediction is not totally unexpected. Since the Kern 
relationship is a more realistic view of the condensate 
action in the tube bundle, one would expect the data to be 
much nearer to the Kern relationship — as it is. The close 
agreement with the Kern relationship also was considered to 
be an indication that the test apparatus were operating 
properly . 

Figure 5.4 shows the normalized, average heat-transfer 
coefficient for up to 30 tubes. The curve generated by these 
data is much smoother than the curve generated by the 
normalized, local heat-transfer coefficients. The value for 
the normalized, average heat-transfer coefficient of the 40th 
tube lies about 26% above the value predicted by Nusselt, 
but lies only 5% below the value predicted by Kern. Assuming 
a generalized form of the average heat-transfer coefficient 
for N tubes, divided by the heat-transfer coefficient for 
the first tube in the bundle, yields: 

h^/h 1 = N _s (5.1) 

Referring to Figure 5.4, the dashed line roughly following 
the data represents a least-squares-f it exponential curve 
for the data, and the exponent derived is s = 0.183. This 
exponent and the exponent suggested by Kern (s = 0.167) are 
in reasonable agreement. 
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Figure 5.4 Inundation Effect on a Smooth Tube Bundle 

(Normalized, Average Heat-Transfer Coefficient) 



2 . With Vapor-Shear Contribution Removed 



Figure 5.3 shows the variation of the normalized, 
local heat-transfer coefficient for up to 30 tubes with the 
contribution due to vapor shear removed. The vapor-shear 
correlation used was the one presented earlier in this report 
(B = 0.834 and n = 0.185). As can be seen from the graph, 
the trend in the data has now become a reversed saw-tooth. 
This is an indication of overcorrecting for the effect of 
vapor shear. This unacceptable trend in the data shows that 
mathematical separation of the vapor-shear contribution is 
not possible as stated by Marto [Ref. 27]. Therefore, 
throughout the remainder of this thesis the vapor-shear con- 
tribution is included, and no further attempt is made to 
remove the contribution due to vapor shear. 

D. THE EFFECT OF INUNDATION ON SMOOTH TUBES WRAPPED WITH 

1.6-MM-O.D. WIRE 

Even though three different pitches (16 mm, 7 . 6 mm and 
4 mm) were tested at this wire diameter, only the graph for 
the optimum pitch will be presented for the normalized, 
local heat-transfer coefficient for a bundle of up to 30 
tubes. The graph for the pitch of 7.6 mm is presented in 
Figure 5.5. The value for the normalized, local heat- 
transfer coefficient for the 30th tube is approximately 170% 
above the value predicted by the Nusselt theory and approxi- 
mately 83% above the value predicted by the Kern 
relationship . 
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Figure 5.5 Effect of 1.6-mm-o.d. Wire on Inundation 

(Normalized, Local Heat-Transfer Coefficient) 



Figure 5.6 shows the normalized, average heat-transf er 
coefficient for a bundle of up to 30 tubes. All three 
pitches are represented, with the values of the normalized, 
average heat-transfer coefficient of the 30th tube for the 
pitches of 16, 7.6 and 4 mm approximately 45%, 118%, and 
120%, respectively, above the value predicted by Nusselt. 

As mentioned previously, the wire wrapping makes two separate 
contributions towards the increase in the heat-transfer 
coefficient — it thins the condensate film and acts as a 
condensate drainage channel. When compared with Figure 5.3, 
Figure 5.5 shows that condensate inundation has less effect 
on wire-wrapped tubes than it does on a set of smooth tubes 
and that a wire pitch of 7.6 mm gives the optimal results. 

The optimum wire pitch or wire diameter is not simply 
based on the largest value of h^/h^ , but also on the actual 
value of h^ . For example, both 4 mm wire pitch and 7.6 mm 
wire pitch resulted in about the same h^/h^ value (see 
Figure 5.6). But the h^ on the tube wrapped with the 7.6 mm 
wire pitch is about 10 percent greater than the tube with 
the 4.0 mm wire pitch. Thus, the 7.6 mm wire pitch gave 
the optimal results for the 1.6-mm-o.d. wire. A more 
thorough discussion of the point is provided later in this 
chapter (Section H) . 

E. THE EFFECT OF INUNDATION ON SMOOTH TUBES WRAPPED WITH 

1.0-MM-O.D. WIRE 

Three different pitches (8 mm, 6 mm and 4 mm) were 
tested at this wire diameter. Again, only the graph for the 
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Figure 5.6 Effect of 1 . 6-mm-o . d . Wire on Inundation 

(Normalized, Average Heat-Transfer Coefficient) 



optimum pitch will be presented for the normalized, local 
heat- transfer coefficient for a bundle of up to 30 tubes. 

The graph for a pitch of 4 mm is presented as Figure 5.7. 

The value for the normalized, local heat- transfer coefficient 
for the 30th tube is approximately 151% above the value 
predicted by the Nusselt theory and 71% above the value 
predicted by the Kern relationship. 

Figure 5.8 shows the normalized, average heat-transfer 
coefficient for a bundle of up to 30 tubes. Again, all 
three pitches are represented. The values of the normalized, 
average heat-transfer coefficient of the 30th tube (for the 
pitches 8, 6 and 4 mm) are approximately 90%, 105% and 117%, 
respectively, above the value predicted by Nusselt. These 
graphs again show that this tube configuration is less 
affected by condensate inundation than the smooth-tube set, 
and that the optimal wire pitch is 4 mm with this smaller 
diameter wire. 

F. THE EFFECT OF INUNDATION ON SMOOTH TUBES WRAPPED WITH 

0.5-MM-O.D. WIRE 

Two pitches (4 mm and 2 mm) were tested at this wire 
diameter. For the normalized, local heat-transfer coeffi- 
cient for a bundle of up to 30 tubes, only the optimum pitch 
of 4 mm will be presented. The graph is presented as Figure 
5.9. The value for the normalized, local heat-transfer 
coefficient for the 30th tube is approximately 101% above 
the value predicted by Nusselt and 37% above the value 
predicted by the Kern relationship. 
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Figure 5.7 Effect of 1.0-mm-o.d. Wire on Inundation 

(Normalized, Local Heat- Transfer Coefficient) 
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Figure 5.8 Effect of 1 . 0-iran-o .d . Wire on Inundation 

(Normalized, Average Heat-Transfer Coefficient) 
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Figure 5.9 Effect of 0 . 5-mm-o .d . Wire on Inundation 

(Normalized, Local Heat-Transfer Coefficient) 



Figure 5.10 shows the normalized, average heat-transfer 
coefficient for a bundle of up to 30 tubes. Both of the 
pitches are represented. The values of the normalized, 
average heat-transf er coefficient for the 30th tube (for the 
pitches 4 and 2 mm) are approximately 71% and 80%, respectively, 
above the value predicted by Nusselt. These graphs also 
indicate that this tube configuration is not affected by 
condensate inundation as much as the smooth-tube set. 

G. THE EFFECT OF INUNDATION ON SMOOTH TUBES WITH A DROPWISE 

COATING 

A set of smooth tubes with a dropwise coating was tested 
primarily to observe the effect of inundation. Figure 5.11 
shows the normalized, average heat-transf er coefficient for 
a bundle of 30 tubes. Also included in the figure for 
comparison are the data for the smooth tube. As is immedi- 
ately evident from the graph, condensate inundation does not 
decrease the performance of the dropwise-coated tubes. The 
ratio of the average outside heat-transfer coefficient of N 
tubes, divided by the heat- transfer coefficient of the first 
tube can be seen to rise at first and then remain steady at 
a value of approximately 1.1. 

Since the coating is non-wetting, no large droplets form; 
rather, the droplets are small when they first form and roll 
down the side of the tube, as soon as they reach a moderate 
size. As these droplets are making their descent down the 
side of the tube, they collect additional droplets until they 
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Figure 5.10 Effect of 0.5-mm-o.d. Wire on Inundation 

(Normalized, Average Heat-Transfer Coefficient) 
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Figure 5.11 Effect of Inundation on a Dropwise Tube Bundle 



reach a point where they are able to fall off the tube. As 
these droplets strike the tubes below, they again collect 
condensate and act as a sweeper of the condensate across the 
surface of the tubes below. 

H . SUMMARY 

Table 3 represents a summary of the results obtained 
from the data sets presented in this thesis. The column 
hi/h Nu is the measure of the local heat- transfer coefficient 
of the first tube compared to the Nusselt prediction. It 
can be considered a measure of the effectiveness of the 
wire-and-pitch combination to thin the condensate film 
between successive wire wraps. Refer to Figure 5.12 for 
clarification of this point. The column h^/h^ is the 
measure of the average heat-transfer coefficient for 30 
tubes compared to the coefficient for the first tube. This 
can be considered to be a measure of the effectiveness of 
the wire-and-pitch combination at handling the condensate 
drainage. The column h^/h^ is the measure of the average 
heat-transfer coefficient for 30 tubes compared to the heat- 
transfer coefficient predicted by the Nusselt theory for a 
single tube. It can be viewed as a measure of the overall 
performance of the 30 tube bundle. The last column (labelled 
s) is the exponent defined in equation (5.1), and is calcu- 
lated using the least-squares fit for each data run. 

As can be seen from the data, it appears that the larger 
the wire, the better it aids in condensate drainage. 
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NOTE: Determined by curve fitting the data for tubes 11 through 30. 




= Exposed surface for diameter 1. 

S 2 = Exposed surface for diameter 2. 

Figure 5.12 Effect of Wire Diameter on Condensation 
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Conversely, the smaller the wire the higher the heat- transfer 
performance of the first tube. These last two observations 
are based on a comparison at a 4 mm pitch but with three 
different wire diameters. The explanation for this is that, 
for a given pitch, as the wire diameter of the wrap is in- 
creased, less and less tube surface area remains exposed to 
the steam flow. In addition to this, the condensate film 
thickness may also increase as the wire diameter is increased. 
Thus, for a given wire pitch, the heat-transfer coefficient 
of the first tube will increase with a decrease in the wire 
diameter. On the other hand, the amount of condensate col- 
lected around the wire increases with increasing wire diameter 
(see Figure 5.12). While the surface-tension forces tend to 
retain liquid at the wire, gravitational forces tend to draw 
liquid along the wire. As the wire diameter increases, the 
gravitational forces increase at a faster rate than the 
surface-tension forces. Therefore, condensate drainage in- 
creases with increasing wire diameter. 

I. OBSERVATIONS 

1 . Smooth Tubes 

During the data runs, filmwise condensation was 
observed with no visible indication of dropwise condensation. 
When the inundation rate was increased, more drops formed 
along the bottom of the tubes , but the droplet size did not 
increase noticeably. However, at the higher inundation 
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conditions there was a noticeable amount of splashing 
occurring in the test condenser. 

2 . Wire-Wrapped Tubes 

As the steam was condensing on the tubes, the conden- 
sate was observed to be drawn to the wire wrap, leaving a 
thin film between the wires (film thickness was dependent 
upon the particular wire diameter and pitch) . The condensate 
then ran down the wire and formed droplets on the bottom of 
the tube. The droplet size was different for each configura- 
tion of wire diameter and pitch. For the larger diameter 
wire, rivulets formed and there was a solid condensate 
"bridge" from the tubes above to the tubes just below them. 

As the inundation rate increased, these bridges became more 
pronounced. It was also noticed that the tendency to splash 
increased as the inundation rate increased, and was quite 
noticeable for all of the tube sets tested. 

3 . Dropwise-Coated Tubes 

As was mentioned previously, the condensate droplets 
from the upper tubes in the bundle act as sweepers of the 
condensate across the surface of the tubes below them in the 
tube bundle. During the data run, the condensate collecting 
on the upper tubes and then falling onto the tubes below was 
sweeping the lower tubes at such a rate that the droplets 
forming on the lower tubes never reached a size larger than 
approximately 2 mm. 
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VI. CONCLUSIONS 



1. The average, outside heat-transfer coefficient for 
30 smooth tubes in a vertical column was 64% of the 
Nusselt coefficient for the first tube in the tube 
bundle. 

2. With wire wrapping, the average, outside heat-transf er 
performance of the tube bundle is considerably improved. 
Significant improvements in the outside heat-transf er 
performance are also possible by utilizing dropwise 
coatings . 

3. The optimum pitch for the 1.6-mm-o.d. wire wrapping 
occurred at 7.6 mm. With this wire arrangement, the 
average, outside heat-transfer coefficient for the 
30 tubes was 92% of the Nusselt coefficient for the 
first tube in the tube bundle. 

4. The best pitch tested for the 1.0-mm-o.d. wire 
wrapping occurred at 4 mm. With this wire arrangement, 
the average, outside heat-transfer coefficient for 

the 30 tubes was 115% of the Nusselt coefficient for 
the first tube in the tube bundle. 

5. The best pitch tested for the 0.5-mm-o.d. wire wrapping 
occurred at 4 mm. In this case, the average, outside 
heat-transfer coefficient for the 30 tubes was 110% 

of the Nusselt coefficient for the first tube in 
the tube bundle. 

6. In simulating a bundle of 30 tubes coated with a 
dropwise promoter, the normalized, average, outside 
heat- transfer coefficient increased by approximately 
10% over the value for the top tube in the bundle. 
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VII. RECOMMENDATIONS 



A. TEST APPARATUS MODIFICATIONS 

The modifications listed below should be considered 
prior to continued use of the test apparatus: 

1. Redesign the test condenser to reduce the possibility 
of the steam bypassing the active test condenser tubes. 

2. Install a larger heating system in the perforated- 
tube water supply tank or modify the existing heating 
system to increase the slow response of the heating 
system. 

The modifications listed below would be beneficial to 
the operation of the test apparatus but are not necessary 
for proper operation: 

1. Redesign the existing controls so that the test 
apparatus is controllable from a single panel. 

2. Redesign the test condenser hotwell to provide a more 
accurate measure of the condensate collection rate. 



B. ADDITIONAL TESTS TO CONDUCT 

The following additional tests would be important in 
the continuation of this investigation: 

1. Conduct tests with additional wire pitches for the 
1 . O-mm-diameter wire to determine the optimum wire 
pitch . 

2. Conduct tests using other dropwise coatings. 

3. Conduct tests with the commercially-available finned 
tubes . 

4. Conduct additional tests to determine the effect of 
vapor shear on the outside heat-transfer coefficient. 
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APPENDIX A 



SAMPLE CALCULATIONS 

A. A sample calculation is performed in this section to 
illustrate the solution procedure used in the data-reduction 
program presented in Appendix C. The calculations are limited 
to the first tube in the bundle only. All thermophysical 
properties were determined from the Tables in Reference 32. 



EXPERIMENTAL CONDITIONS : 

Pressure Condition 

Inundation Condition 

Inlet Temperature of Cooling Water 

Outlet Temperature of Cooling Water 

Saturation Temperature 

Cooling Water Rotameter Setting 



Atmospheric 
5 tubes 
24.10 °C 
28.32 °C 
100.65 °C 
21.7% 



1 . Determination of Average Bulk Temperature 



= (T 



ci 



T )/2 
co 



T = (24.10 + 28 . 32) /2 

b 



T, = 26.21 °C 

b 

2. Thermophysical Properties (evaluated at T^) 

P = 5.83 

r 

P c = 997 kg/m 3 
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U c = 855 xio 6 N-s/m 2 



C pc = 4.179 kJ/kg-K 



-3 

k c = 618 x 10 W/m-K 



m = 0.243 kg/s 

c 



4 . Determination of Cooling Water Velocity 



V = m / ( p -A. ) 
c c c 1 



V = (0.243)/ (997 ) • (1.56 x 10 4 ) 

c 



V c = 1.56 m/s 



5 . Determination of Reynolds Number 



Re = (p c V c D.)/y c 



Re = (997) (1.56) ( .0141)/855 x 10 



-6 



Re = 25,649 



6 . Determination of Heat Transfer 

Q = m (T - T . ) C 

c co ci pc 

Q = 0.243(28.32 - 24.1)4179 

Q = 4285 W 
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7. Determination of Heat Flux 



q = 


Q/(tt D q L) 


q = 


4285/tt (0.015875) (0.305) 


q = 


281,700 W/m 2 



8. Determination of Nusselt Coefficient 



h 

Nu 


"k 3 p 2 h. (9.81) 1 1/3 

= 0.651 — 

y^ D q 
L f o ^ 



Assume : 





T = T 

f sat 

T f = 100.65 °C 


p f 


= 957.8 kg/m 3 


h fq 


= 2255 J/kg 


k f 


= 680.52 W/m-K 


P f 


= 277.6 x 10" 6 N-s/m 2 


h 

Nu 


„ r (680.52xl0 -3 ) 3 (957.8) 2 (2255xl0 3 ) (9.81), 1 / 3 

— U . DD1 l ~ J 

(277.6x10 ) (0.015875) (281,700) 


h 

Nu 


= 11,243.6 W/m 2 • K 
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9. Determination of T 



f ,c 



T j- = T . - q/ ( 2 • h ) 

f/C sat Nu 



T = 100.65 - 281,700/2(11,243.6) 

r , c 



T = 88.12 °C 

f ,c 



10 . Thermophysical Properties 



k f = 674.6 x 10 3 W/m-K 



= 966 kg/m' 



= 320.3 x 10 ^ N*s/m 2 



h fg = 2288.3 x 10 J J/kg 



11. Determination of Nusselt Coefficient 



Nu 



= 0.651 



k f 4 h fq 9 - 81 

p £ D o 5 



1/3 



Nu 



0 r (674.6x10 3 ) 3 (966) 2 (2288.3xl0 3 ) (9.81) , 1/3 
0.651 l — J 

(320.3xl0 _b ) (0.015875) (281,700) 



Nu 



= 10,739.6 W/m K 



12. Determination of T 



f,c 



T r = T . - q/(2*h„ ) 

f , c sat Nu 
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100.65 - 281,700/(2) (10 ,739 .6) 



f ,c 



T f = 87.53 °C 

f ,c 



13 . Determination of Logarithmic Mean Temperature 



T - T . 

LMTD = (T - T .)/£n(=^ 

co ci T ^ - T 

sat co 



LMTD = (28.32 - 24 . 10 ) An ( ; ff 1^° ) 



LMTD = 74.42 °C 



14 . Determination of Overall Heat-Transfer Coefficient 

U Q = q/LMTD 

U q = 281,700/74.42 

U q = 3785.3 W/m 2 • K 

15. Determi nation of Inside Heat-Transfer Coefficient 
Assume : 

h. 



C f = 1.1 



C. = 0.028 

l 



C. Re 0 * 8 p r ° * 33 3 Cf 
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h i = 6l4l (° • 028 ) (25, 649) °- 8 (5.83) 0,333 (1.1) 
hu = 8181 W/m 2 K 

16 . Determination of Inner Wall Temperature 

T = T, + q D /(h. D. ) 

w b ^ o l l 

T = 26.21 + (281,700) (0.015875)/ (8181) (0.0141) 

w 

T = 64.98 °C 



17 . Determination of y. ^ at the Average Wall Temperature 

y = 433 x 10 ^ N*s/m 2 

w ' 



18. Determination of Correction Factor 



C- = (y /y ) 
f ,c c' w 



0.14 



C f,c = 1 * 1 



19. Determination of Outisde Heat-Transfer Coefficient 



_±_ _ ^ _ R 

U D. h. w 
o 11 



0.015875 



3785.3 (0.0141) (8181) 



- 0.000042925 



h = 11,957 W/m K 

o 
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APPENDIX B 



UNCERTAINTY ANALYSIS 



The uncertainty analysis performed for this thesis is 
based on the development described in Kline and McClintock 
[Ref. 33] . Kanakis [Ref. 12] employed this development and 
derived the expressions for uncertainty in the variables of 
this thesis also. In this thesis, the results presented by 
Kanakis [Ref. 12] in Appendix C will be utilized, but where 
necessary, corrections will be pointed out. 



A. UNCERTAINTY IN THE COOLING WATER VELOCITY 



6 V 

c 

V 

c 



6m 6p 

(— ) + (— 



m 



„ 6A. 0 1 1/2 

V + HrV 



A. 

1 



with the following uncertainties assigned: 



6m = 

c 


+ 


0.01 kg/s 


II 

<J 

CL 


+ 


3 kg/m^ 


6A. 


+ 


0.00000022 m 



B. UNCERTAINTY IN THE REYNOLDS NUMBER 



6 Re 
Re 



^P c 2 

(-r 2 -) + ( 



6V c 2 6D i 2 6y c 2 1 1/2 
C ) + (-FT^r + (^r 2 -) 



V 



D. 

l 



y. 
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with the following uncertainties assigned: 



<Sp c = ± 3 k 9/ m ' 



6 V as calculated in Section A 
c 



6D. = ± 0.0001 m 

l 



6y = ± 10 xio 6 N-s/m 2 

c 



C. UNCERTAINTY IN HEAT TRANSFER 



6Q = 
Q 



6m 0 6T „ 6 T . ~ 

(— ) + ( n T— S 9 --) + (m ^-) + ( 



m 



T -T . 
co ci 



T -T 
CO ci 



6C 



^2 

'pc 



1/2 



with the following uncertainties assigned: 



6m c = ±0.01 kg/s 



6T = ± 0.025 °C 

co 



6T . = ± 0.025 °C 

ci 



SC = ± 2 J/kg-C 

pc 



D. UNCERTAINTY IN THE HEAT FLUX 



= 

q 



( ^ o )2 + ( “^ )2 + ( ^ r )2 



1/2 



with the following uncertainties assigned: 
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6Q as calculated in Section C 

<SD = ± 0.0001 m 

o 

6L = ± 0.001 m 



E. UNCERTAINTY IN h XT 

Nu 



6h 



Nu 

*Nu 



6k 



f ? 2 

— — ) + ( — 

k ' 3 



6p 



f -) 2 



i <5h o 

+ ( T TT7 2 ’ 2 
fg 



(i M ) 2 

3 q 1 



.1 6 y f >2 1 6 D o ,2 

+ ( 3 U7 ) + ( 3 ~D + 



/I 5q^ 2 
l 3 



1/2 



with the following uncertainties assigned: 
5k f = ± 0.0005 W/m-K 

6p f = ±3 kg/m^ 



6h 



fg 

6g 






5D 

o 

6q 



= ±0.12 J/kg 

= ± 0.0005 m/s 

— 6 2 

= ± 4.9 x 10 N's/m 

= ± 0 . 0001 m 

as calculated in Section D 
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F. UNCERTAINTY IN T 



f ,C 



ST 



f,c 



■f ,C 



ST „ . 0 Sh„ „ 

/ sat v 2 ,6q^2 . / NU) 2 

' T _ . ' K q 1 1 h, ’ 



sat 



Nu 



1/2 



with the following uncertainties assigned: 



ST . = ± 0.025 °C 

sat 



Sq as calculated in Section D 



5h„ as calculated in Section E 
Nu 



G. UNCERTAINTY IN LOGARITHMIC-MEAN-TEMPERATURE DIFFERENCE 



SLMTD 

LMTD 



ST (T -T .) 
sat co ci 



(T sat- T ci )(T sat- T co )tn( T Sat -T Cl 



sat co 



ST 



co 



(T sat- T co Un( T Sat -T Ci ) 



sat co 



ST 



ci 



T -T . 

(T -T ) Jin ( sa ^ £A) 

U sat i ci ,x 'T -T 1 



sat co 



1/2 



with the following uncertainties assigned: 
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<5T = 

sat 


± 0.025 °C 


ST 

CO 


± 0.025 °C 


ST . 

Cl 


± 0.025 °C 



H. UNCERTAINTY IN OVERALL HEAT-TRANSFER COEFFICIENT 



6U 

o 

U o 


, Sq. 2 , , 5LMTD. 2 1 

( q ) + ( LMTD ) 


with the following 


uncertainties assigned: 



5q as calculated in Section D 

5LMTD as calculated in Section G 



I . UNCERTAINTY IN 


INSIDE HEAT-TRANSFER COEFFICIENT 


II 

<o 


, 5k c,2 , ^ D i . 2 0.8 SRe.2 ,0.333 SPr . 2 

1 k ; 1 D. 1 Re ; 1 Rr ' 

C 1 




SC. 2 0.145 (u /u ) .I 1/2 

+ ( + ( , C W ) 2 
C . y /y 

1 H c p w J 



with the following uncertainties assigned: 



+1 

II 

a 

to 


0.0005 W/m*K 


+1 

II 

•H 

Q 

to 


0.0001 m 
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SRe as calculated in Section B 



If 

u 


± 0.063 


6C. = 

l 


± 0.0001 




-a. o 



S(y /y ) = ±4.9x10 N*s/m 

c w 



J. UNCERTAINTY 


IN TEMPERATURE DIFFERENCE 


6DT 

DT 


x ~ 6h. „ SD ~ 

(^) + (-r^) 2 + (-R^) 2 

_ q h. d o 



with the following uncertainties assigned: 



5q as 


calculated in Section D 


5h . as 

l 


calculated in Section I 


SD 

o 


± 0.0001 m 


6D. 


± 0.0001 m 



K 



UNCERTAINTY IN OUTSIDE HEAT-TRANSFER COEFFICIENT 




with the following uncertainties assigned: 

6U as calculated in Section H 

0 

5R = ± 0.00001 m 2 • K/W 

w 

6h. as calculated in Section I 

1 

For the data presented in Appendix A, the following uncertain- 
ties are the result: 



quantity 


uncertainty 


V = 1.56 

c 


± 0.064 m/s 


Re = 25,649 


± 1x12 
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quantity 



Q 


— 


4258 


q 


= 


281,700 


h. T 

Nu 


= 


11,243.6 


f ,c 


= 


88.12 


LMTD 


= 


74.42 


U 

O 


= 


3785.3 


h. 

l 


= 


8181 


DT 


= 


38.77 


h 

o 




11,957 



uncertainty 
± 169 w 

± 12,010 W/m 2 

± 163.4 W/m 2 • K 

± 4.0 °C 

± 0.625 °C 

2 

± 164.48 W/m -K 

± 292.6 W/m 2 • K 

± 2.19 °C 

± 2287 W/m 2 • K 
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APPENDIX C 



DATA REDUCTION AND PLOTTING PROGRAMS 



1000 ! 
1005 ! 
1010 ! 
1015 
1020 
1025 
1030 
1035 
1 040 
1045 
1050 
1055 
1060 
1065 
1070 
1075 
1080 
1085 
1090 
1095 
1100 
1105 
1110 ! 
1 1 15 ! 
1120 ! 
1125 ! 
1130 
1135 
1 140 
1145 ! 
1 150 ! 
1155 ! 
1160 ! 
1165 
1170 
1175 
1180 ! 
1185 ! 
1190 ! 
1 1 95 
1200 ! 
1205 ! 
1210 
121 5 
1220 
1225 
1230 
1235 ! 
1240 ! 
1245 
1250 ! 
1255 ! 
1260 
1265 
1270 ! 
1275 
1280 
1285 
1290 



FILE NAME : DRP3 
DATE: May 9. 1984 

BEEF’ 

PRINTER IS 1 

PRINT USING "4X . * ,M Se J ec t Option: 

PRINT USING "6X, ,,M 0 Taking data or re-processing previous data 

PRINT USINC "6X . **** 1 Plotting previous data 

PRINT USING ”6X.*'**2 Computing exponent for experimental data 

PRINT USING **6X ♦ ”*'3 Plottinq on LOG-LOG 

PRINT USING ”6X. #,M 4 Labelling 

INPUT lop 

IF lop=0 THEN CALL Main 

IF lop* 1 THEN CALL Plot 

IF Iop=2 THEN CALL Expo 

IF Iop=3 THEN CALL Lp lot 

IF Iop*4 THEN CALL Label 

END 

SUB Main 
COM /Ci/ C< 7 ) 

DIM Tci<2) ,Tco<4,2) . Ti<4) .Mft(4) .Vu<4> .Ho<4) 

DIM To(4).Ts<1 ) ,Tb(4) ,R3(4).R4<4) , S3<4) ,S4<4) ,Fc<4) ,Mfi(4).Mfs<4> 

ASSIGN COEFFICIENTS FOR THE 8-TH ORDER 
POLYNOMIAL FOR TYPE-T ( COPPER-CONSTANTAN ) 

THERMOCOUPLES 

DATA 0.10086091 .25727.94369.-767345.8295.78025595.81 
DATA -9247486589. 6. 97688E+1 1 ,-2.66192E+13,3.94078E+1 4 
READ C(*> 

ASSIGN CONSTANTS FOR ROTAMETER CALIBRATION 
LINES (kg/m in) 

DATA -1 .256E-4.5.866E-3.2.456E-4.-3.352E-3.-5.878E-4 
DATA 1 .121E-2.1 .237E-2.1 .221E-2.1 .233E-2 J .232E-2 
READ Mf i < * ) . M f s( * ) 

ASSIGN SIEDER-TATE COEFFICIENT AND EXPONENT 
FOR REYNOLDS NUMBER 
Ex= . 8 



ASSIGN GEOMETRIC VARIABLES 
Di=.0141 ! Inner diameter <m) 

Do=. 015875 ! Outer diameter <m) 

Ktm=21 .9 ! Thermal conductivity of titanium 

L*.305 ! Condensing length <m) 

Pt = Do*1 .5 ! Transverse tube p i tch- to-d lairie ter 



(W/m-K) 
rat i o 



COMPUTE THE MEAN STEAM FLOW AREA IN THE TEST CONDENSER <m 2) 
Af = (Pt* 2-PI*Do 2/4 ) /Pt«2*L 

COMPUTE INSIDE AREA AND WALL RESISTANCE 
Ai-PI*Di ' 2/4 

Rui = Do*L0G( Do/Di >/<2«Ktm> 



PRINTER IS 70i 
CLEAR 70S 
BEEP 

INPUT “ENTER MONTH. DATE. AND TIME < MM : DD : HH : MM : SS >*’ . T i me$ 
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1295 

1300 

1305 

1310 

1315 

1320 

1325 

1330 

1335 

1340 

1345 

1350 

1355 

1360 

1365 

1370 

1375 

1380 

1385 

1390 

1395 

1400 

1405 

1410 

1415 

1420 

1425 

1430 

1435 

1440 

1445 

1450 

1455 

1460 

1465 

1470 

1475 

1480 

1485 

1490 

1495 

1500 

1505 

1510 

1515 

1520 

1525 

1530 

1535 

1540 

1545 

1550 

1555 

1560 

1565 

1570 

1575 

1580 

1585 

1590 

1595 



OUTPUT 709 : "TD" : 1 ime$> 

Sgy-0 
Sgf *0 
Sgyf *0 
Sgf 2*0 
BEEP 

INPUT "VAPOR SHEAR DATA ( 1 - Y . 0-N ) ?" . Ivd 

IF Ivd'O THEN 

BEEP 

INPUT "ENTER OUTPUT MODE ( 1 -SHORT , O-LONG >". Jop 

END IF 

BEEP 

INPUT "ENTER THE INPUT MODE < 1 -3054A .2-FILE >". Im 

IF Im-2 THEN 

BEEP 

INPUT "ENTER THE NAME OF THE EXISTING DATA F ILE" . OlddataS 

PRINT USING "10X,""This analysis is for data file 1 OA" ; OlddataS 

BEEP 

INPUT "ENTER 1 IF PROCESSING 1983 DATA" . Idt 
ASSIGN ©File2 TO OlddataS 
END IF 

IF Im-1 THEN 
BEEP 

INPUT "GIVE A NAME FOR THE DATA FILE TO BE CREATED" . NewdataS 

CREATE BDAT NeudataS,25 

ASSIGN ©Filel TO NeudataS 

END IF 

BEEP 

INPUT "GIVE A NAME FOR THE OUTPUT FILE" ,F i le_out$ 

BEEP 

INPUT "ENTER TUBE TYPE < 0 = PLAIN . 1 = ROPED)" . I tt 
IF Itt-0 THEN 

PRINT USING " 1 OX . ""Tube type : Plain 

Ci«. 028 
END IF 

IF Itt=1 THEN 

PRINT USING " 1 OX , ""Tube type : Roped 

Ci*. 057 
END IF 

PRINT USING " 1 OX . ""S ieder-Tate constant - "" . Z . 4 D" : C i 
BEEP 

INPUT "ENTER THE INUNDATION CONDITION <1=5 TUBES. 2=30 TUBES)". Mi 
IF Mi=2 THEN PRINT " Inundation condition : 30 TUBES" 

IF M i = 1 THEN PRINT " Inundation condition : 5 TUBES" 

IF I v d = 0 THEN 
BEEP 

INPUT ’’WANT TO INCLUDE VAPOR SHEAR < 1 = Y . 0 = N> 7 " . Ivs 

IF Ivs-1 THEN 

BEEP 

INPUT "ENTER B-VALUE FOR VAPOR SHEAR CORRELATION" . B 
BEEP 

INPUT "ENTER EXPONENT FOR VAPOR-SHEAR CORRELAT ION" . Nvs 
END IF . 

END IF 

CREATE BDAT File_outS.6 
ASSIGN ®File3 TO File_outS 
Ja=0 
Nrun=0 

FOR 1=0 TO 4 
S3< I >=0 , 

S4< I > = 0 . 
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1600 NEXT I 
1605 IF Im-1 THEN 
1610 BEEP 

1615 INPUT "ENTER MANOMETER READINGS (HLW ,HLM , HRW , HRM)" .Hlu.Hlm.Hru.Hrm 
1620 OUTPUT SFilel iHlw.Hlm.Hru.Hrm 
1625 ELSE 

1630 IF Idt =0 THEN ENTER 9F 1 1 e2 : HI u . H lm . Hru .Hrm 
1635 END IF 

1640 Dp = Hlm-Hrm + (Hlui-Hlm-Hrui+Hrm)/13.5 
1645 Mdot *FNMdot (Dp) 

1650 Nrun'Nrun* 1 
1655 OUTPUT 709 : "TD" 

1660 ENTER 709:TimeS 

1665 IF lvd = 0 AND Jop=0 THEN PRINT •• •• 

1670 IF Nrun- 1 THEN 

1675 PRINT USING " 1 OX . ""Month . date, and time: . 1 5A" : T imeS 

1680 PRINT 

1685 END IF 

1690 IF Im = 2 THEN Rdf 

1695 BEEP 

1700 INPUT "ENTER FLOW METER READINGS (AS PERCENTAGES)” ,Fm1 ,Fm2 
1705 DISP "START COLLECTING CONDENSATE" 

1710 BEEP 
1715 WAIT 20 

1720 OUTPUT 709 : ”AR AF0 AL19" 

1 725 OUTPUT 722 : "F 1 R1 T1 Zl FL1 " 

1 730! 

1735! READ INLET WATER TEMPERATURES 
1 740! 

1745 FOR 1=0 TO 2 
1750 OUTPUT 709 ; "AS SA" 

1755 ENTER 722:Tci<I) 

1760 CALL Tvsv( Tc i ( I ) > 

1765 Tc i ( I ) =FNTemp < Tc i < I ) . I ) 

1770 NEXT I 
1 775 ! 

1780! READ OUTLET WATER TEMPERATURES 
1785! 

1790 1 1-2 

1795 FOR 1=0 TO * 

1800 IF 1=0 OR 1=3 THEN 
1805 In- 2 

1810 else 

1615 Iu=1 

1820 END IF 

1825 FOR J = 0 TO Iu. 

1830 Ii-Iid 

1835 OUTPUT 709 ; ”AS SA" 

1840 ENTER 722:Tco<I.J) 

1845 CALL Tvsv< Tco( I , J) > 

1850 Tco(I ,J)=FNTemp(Tco<I,J> ,Ii) 

1855 NEXT J 
1860 NEXT I 
1 00S * 

1870! READ STEAM TEMPERATURES 
1875! 

1880 FOR 1=15 TO 16 
1885 OUTPUT 709 : "AS SA" 

1890 ENTER 722:Ts<I-15) 

1895 CALL Tvsv<Ts<I-15)) 

1900 Ts(I- !5)=FNTemp<Tsd-15) .I> 
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1905 NEXT I 
1910! 

1915! READ CONDENSATE TEMPERATURE 
1920! 

1925 OUTPUT 709 ; "AS SA" 

1930 ENTER 722:Tcon 
1935 CALL Tvsv(Tcon) 

1940 Tcon-FNTemp< Tcon , 1 7 > 

1945! 

1950! READ VAPOR TEMPERATURE 
1955! 

1960 OUTPUT 709 ; "AS SA" 

1965 ENTER 722;Tv 
1970 CALL Tvsv(Tv) 

1975 Tv=FNTemp( Tv .18) 

1980! 

1985! READ VAPOR PRESSURE 
1990! 

1995 OUTPUT 709 : "AS SA" 

2000 ENTER 722:P_volts 
2005! 

2010! COMPUTE AVERAGE WATER TEMPERATURES AT INLFT 
2015! 

2020 Ti<0>=Tci<0> 

2025 T i < 1 > = <Tci<0) + Tcj< 1 > >*.5 
2030 Ti(2>=Tci<1) 

2035 Ti<3)=(Tci(1>+Tc]<2>>«.5 

2040 Ti(4)=Tci(2> 

2045! 

2050! COMPUTE AVERAGE WATER TEMPERATURES AT OUTLET 
2055! 

2060 FOR 1 = 0 tq zj 
2065 IF 1=0 OR 1=3 THEN 

2070 To<I>-<Tco<I.0>*Tco<I.1 >+Tco< I .2) >* . 3333 

2075 ELSE 

2080 Tq<I>=(Tco<I.O)+Tco(I. 1 >>*.5 
2085 END IF 
2090 NEXT I 

2095 Tsa=<Ts<OMTs< 1 >>* .5 
2100 Pvap=FNPvsv < p_vo 1 t s ) 

2105 Tsat =FNT vsp< Pvap* 1 33 . 322 ) 

2110 Dsup=Tv-Tsat 
21 15! 

2 1 20» READ INFORMATION FOR CONDENSATE FLOW RATE 
2*1 25 ! 

2130 IF lvd = 0 THEN 
2135 BEEP 

2140 INPU T "ENTER INITIAL AND FINAL LEVELS IN HOT WELL T\H1.H2 
2145 Bh=H2-H1 

2150 IF Nrur, MOD 5=1 THEN Msun. = 0 
2155 Mf 1=540. 4836*Dh 

2 * 6G Mdl =Mf 1 *FNRhou< Tsat- 1 0)* 1 .0E-6/60 
2165 Msur> = Msain+Mf 1 

2170 IF Mi =2 AND Nrun<>30 AND Nrur. MOD 5 = 0 THEN 
2175 Mave = Msurt/5 

2180 Set® (Mave*KNRhow< Tsat- 1 0>/1 0 6+ . 03238 > / . 042 1 32 

2185 END IF 

2190! 

2195 Rdf: ! 

2200 ! 

2205 IF lvd=0 AND Jop = 0 THEN 



91 



2210 PRINT USING "10X.""Run number * .DD" :Nrun 

2215 PRINT " Tube » : 1 2 3 4 5' 

2220 END IF 
2225 END IF 

2230 IF Ivd* 1 AND Nrun-1 THEN 

2235 PRINT USING ” 1 0X .""Data * Vv Retp Nuc Y F 

2240 END IF 
2245 IF Itn-2 THEN 

2250 IF Nrun MOD 5-1 AND Mi-2 AND Nrun>5 THEN ENTER ®File2:Fpt 
2255 ENTER $F i le2 : T i (• > , To( * > . Tsa, Tcon . Tv .Pvao , Tsat , Dsup ,Fm1 ,Fm2 
2280 ENTER 9File2:H1 ,H2 
2265 END IE 

2270 IE lvd=0 AND Jop-0 THEN 

2275 PRINT USING "1 OX .""Inlet temp (Deg C> : .5<DDD . DD.2X)" : Ti (« > 

2280 PRINT USING "1 OX .""Out let temp (Deg C> .5(DDD .DD.2X >" ; To(« > 

2285 PRINT USING " 1 OX , ""Sat urat ion temperature - "".3D.DD."" (Deg C> :Tsat 

2290 PRINT USING "10X. ""Degree of superheat - "".3D.DD,"" (Deg C) ;Dsup 

2295 PRINT USING "1 OX . ""Stat ic pressure - "" . 3D .DD ."" (mm Hq) :Pvap 

2300 END IF 
2305! 

2310! CALCULATE AVERAGE BULK TEMPERATURES 
2315! 

2320 IF lvd=0 THEN Nx = 4 

2325 IF Ivd- 1 THEN Nx = 0 

2330 FOR 1=0 TO Nx 

2335 Tb(I)=(Ti(I)+To(I)>».5 

2340 NEXT I 

2345! 

2350 IF M i = 1 OR (Mi-2 AND Nrun<6> THEN As1=0. 

2355 IF M i=2 AND Nrun>5 THEN S'=As1 

2360 SI =As 1 

2365 FOR J=0 TO Nx 

2370 IF J=0 THEN Cwf=Fm1 

2375 IF J=1 THEN Cuf=Fm2 

2380 Mf*Mf i ( J > +M f s( J > *Cwf 

2385 T x = T b ( J > 

2390 Vu(J)=Mf/(FNRhow(Tx)«Ai ) 

2395' 

2400! CALCULATE INSIDE AND OUTSIDE COEFFICIENTS 
2405! 

2410 Rew=FNRhoiii< T x ) * Vu< J ) *D i /FNMuui ( T x ) 

2415 Cf=1. 

2420 Q-Mf*FNCpui(Tx)*(To(J)-Ti(J)> 

2425 IF Nrun=1 AND J=0 THEN 01=0 
2430 Qp=Q/(PI»Do«L> 

2435 IF ( Mi = 1 OR (Mi=2 AND Nrun<6)) AND J=0 THEN 
2440 Tf i lm=Tsat 
2445 K.f = FNKu< Tf i lm> 

2450 Rhof =FNRhou( Tf ilm) 

2455 Hfg=FNHfg(Tsat>«1000 
2460 Muf = FNMuu(Tf i lit.) 

2465 Hnu= .65 1 *Kf * < Rhof ' 2*Hf g*9 . 799/ (Muf *Do*0p) > .3333 
2470 Tf i Imc-Tsat-Op/Hnu* .5 

2475 IF ABS((Tf ilmc-Tf ilm)/Tf ilmc)>.01 THEN 
2480 T f i lm= T f i lmc 
2485 GOTO 2445 
2490 END IF 

2495 IF I vd=0 AND Jop=0 THEN 

2500 PRINT USING "1 OX . ""Mussel t coefficient for first tube = "".5D.D."" (W/m 2 

K ) :Hr,u 

2505 END IF 
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2510 END IF 

2515 IF J-0 AND lvd=0 AND Jop-0 THEN 
2520 IF Mi- 1 AND Nrun-1 THEN Hol-0. 

2525 PRINT " Tube V u Heat flux Cona coef R1 R2 

RR*’ 

2530 PRINT *• xr fm/s> <W/m 2) (W/m 2.K>" 

2535 END IF 

2540 IF Nrun-1 AND Ivs=1 AND lvd*0 THEN 
2545 Vg*FNVvs t < Tsat ) 

2550 Vv=Mdo t * Vg/Af 
2555 Dt-<Tsat‘Tf i lm>*2 
2560 Muf-FNhuu<Tf ilm) 

2565 Hfg=FNHfg< Tsat >*1000 
2570 Kf-FNKu<Tf ilm) 

2575 Rhof«FNRhou<Tf ilm> 

2580 Re tp“Rhof *Vv*Do/Muf 

2585 Ff-9.799*Do*Muf*Hfg/<Vv'2*Kf*Dt> 

2590 IF I vd=0 THEN 

2595 Nuc=B*Ff ' Nvs*Retp .5 

2600 Hoc*Nuc*Kf /Do 

2605 Tf i lmc=Tsat -Qp/Hoc* . 5 

2610 IF ABSCTf i lm-Tf i lmc >> . 1 THEN 

2615 Tf i lm-Tf i lmc 

2620 GOTO 2555 

2625 END IF 

2630 END IF 

2635 END IF 

2640 Muw = FNf1utii ( T x ) 

2645 Hi=FNKu(T x )/Di*Ci*Reui Ex* ( FNPru < Tx ) > ‘ .3333*C* 

2650 Dt=Op/H i*Do/D i 

2 655 Cfc = (f1uu/(FNMuw<Tx*Dt) ) ) “ . 14 

2660 IF ABS< <Cf-Cfc>/Cfc>> .01 THEN 

2665 Cf=<Cf +Cfc>*.5 

2670 GOTO 2645 

2675 END IF 

2680 Lmtd=<To<J>-Ti<J)>/LOG<<Tsat-Ti<J))/<Tsat-To<J)>> 

2685 IF Nr un = 1 AND Ivs-1 AND lvd = 0 THEN 

2690 Fc< J) =Hoc/Hnu* < Q/O 1 ) .3333 

2695 Mdot “Mdot -Q/Hf g 

2700 END IF 

2705 Uo=Qp/Lmtd 

2710 Ho<J)=1/( 1/Uo-Do/(Di*Hi)-Ru) 

2715 IF I vd= 1 THEN 
2720 Dt=0p/Ho<0) 

2725 Tfilm=Tsat-Dt/2 
2730 Vg=FNVvs t < Tsat ) 

2735 Rhof =FNRhou< T f i lm ) 

2740 Vv=Mdot*Vg/Af 
2745 Muf = FNMuu< T film) 

2750 Kf = FNKw< T f 1 1 m ) 

2755 Hfg=FNHfg< Tsat >*1000 

2760 Re tp-Rhof *Vv*Do/Muf 

2765 Ff=9.799-Do*Muf*Hfg/<Vv' 2*Kf*Dt> 

2770 Nuc-Ho(0)*Do/Kf 
2775 Y=LOG<Nuc/Retp' .5) 

2780 F =L0G< Ff > 

2785 Sgy=Sgy 4 Y 
2790 Sqf -Sgf + F 
2795 Sgyf-SgyF+Y*F 
2600 $ Q f 2=Sqf 2 + F 2 
2805 END IF 
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2b1U if lva-0 THEN 

2815 IF Ivs-1 THEN Ho( J> «Ho< J> /Fc < J ) 

2820 Rr-Uo/Ho(J) 

2825 SI-SUHo(J) 

2830 IF Nr an HOD 5=1 THEN 

2835 IF hi-1 OR (Hi =2 AND Nran=31 > THEN Ja*0 
2840 IF Hi=2 AND 5<Nrun AND Nran<30 THEN Ja=Nrun-1 
2845 IF Mi=2 AND 35<Nran THEN Ja-Nran-1 
2850 END IF 

2855 IF Hi - 1 OR <30<Nrun AND Nran<36 AND Hi - 1 > OR Nran<6 THEN 
2860 Rl-Ho( J)/Ho(0> 

2865 R2-S1 /< ( J+1+Ja>*Ho<0> ) 

2870 ELSE 

2875 R1-Ho(J)/Ho1 

2880 R2-S1 /( < J+1*Ja>*Hol > 

2885 END IF 
2890 r 

2895! PRINT RESULTS 
2900! 

2905 IF Jop = 0 THEN 

2910 PRINT USING " 1 1 X . DD . 4X , DD . DD . 2X , 2( D . 5DE >?X > . 3( Z . 4D , 2X > M ; J+ 1 ♦ Ja .Vw< J ) . Q P . Ho 
< J> ,R1 , R2 . Rr 

2915 IF Im-1 AND J=4 THEN 
2920 BEEP 

2925 INPUT "OK TO ACCEPT THIS DATA SET < 1 «Y . 0«N>r\Qlcs 

2930 IF Of s<> 1 THEN 

2935 Nran=Nran-1 

2940 GOTO 1650 

2945 END IF 

2950 END IF 

2955 END IF 

2960 ! 

2965 IF Hi =2 AND Nran<6 AND J=0 THEN 

2970 Ho1=Ho1+Ho<0>/5 

2975 END IF 

2980 FOR K = 0 TCi 4 

2985 IF K = J THEN S3<K >=S3<K ) + Rl 

2990 IF K= J THEN S4 <K> =S4(K > +R2 

2995 NEXT K 

3000 END IF 

3005 NEXT J 

3010 IF lvd = 0 THEN 

3015 IF Nr an HOD 5=0 THEN 

3020 FOR K=0 TO 4 

3025 R3(K > = S3< K > /5 

3030 R4(K)=S4(K)/5 

3035 S3<K ) =0 . 

3040 S4<K > =0 . 

3045 NEXT K 

3050 IF Hj =2 AND Nran HOD 5=0 AND Nran<>30 THEN As 1 =Nran*R4 < 4 >«Ho 1 

3055 1 

3060! PRINT AVERAGE RATIOS 
3065* 

3070 IF Jop=0 OR <Jop=1 AND Nran<6> THEN 
3075 PRINT 

3080 PRINT " Tube * R3 R4“ 

3085 END IF 
3D90 r GR J-1 TO 5 

3095 PRINT USING ’* 1 2X , DD . 2< 4X . Z . 4D > M ; J+ Ja . R3< J- 1 > , R4 < J- 1 ) 

3100 OUTPUT $Fi le3:J+Ja.R3<J-1 >.R4(J-1 ) 

3 1 0 5 NEXT J 
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3110 END IF 

3115 IF N run MUD 5*0 AND Mi-? AND Nrun<>30 AND Im-1 THEN 
3120 BEEP 

3125 PRINT USING "10X.""Set porous-tube floumeter reading to "" t 3D . D , "" PERCENT 

: Set 

3130 END IF 
3135 END IF 
3 1 40 IF Ivd-l THEN 
31 45 Ey*EXP(Y> 

3150 PRINT USING "1 2X . DD . 5X . Z . DD .2X , Z . 3DE .2X , 3D . DD . 2X . 2< Z . 3D .2X >" ;Nrun , Vv ,Retp , 
Nuc , Ey . F f 

3155 IF Im-1 THEN 
3160 BEEP 

3165 INPUT "OK TO ACCEPT THIS DATA SET < 1 - Y , 0-N > 7" ,0k s 

3170 IF OksOI THEN 

3 1 75 Nrun = Nrun- 1 

3180 GOTO 1605 

3185 END IF 

3190 END IF 

3195 OUTPUT $Fi le3;Ff ,Ey 

3200 END IF 

3205! 

3210 IF Nr un MOD 5-1 AND Mi-2 AND Nrun>5 AND Im-1 THEN 

3215 OUTPUT ©File! :Set 

3220 Mpt=-8. 36 1613+10. 076742* Set 

3225 END IF 

3230 IF Im-1 THEN 

3235 OUTPUT ®F i le 1 ;Ti(*),To<*).Tsa.Tcon,Tv. Pvap , Tsat . Dsup .Fml , Fm2 
3240 OUTPUT GFilel : HI ,H2 
3245 BEEP 

3250 INPUT "WILL THERE BE ANOTHER RUN ( 1 - Y . 0 = N > 7" , Go_on 

3255 IF Go_on= 1 THEN 

3260 IF lvd=0 THEN 1650 

3265 IF Ivd-1 THEN 1605 

3270 END IF 

3275 ELSE 

3280 IF Mi =2 AND Nrun<30 THEN 
3285 IF I vd=0 THEN 1650 
3290 IF Ivd-1 THEN 1605 
3295 END IF 

3300 IF Mi* 1 AND Nrun<10 THEN 
3305 IF I vd=0 THEN 1650 
3310 IF Ivd-1 THEN 1605 
3315 END IF 
3320 END IF 

3325 IF Im- 1 THEN PRINT USING "10X.DD."” Data runs were stored in file . 1 0A" ; 
Nr un , NeudataS 

3330 PRINT 

3335 PRINT USING M 1 OX . M "P Jot data are stored in file "" . 1 4A" : F i le_out$ 

3340 IF Ivd-1 THEN 

3345 Nvs= < Nr un*$gy f -Sgy*Sgf ) / < Nrun*Sgf 2-Sgf 2> 

3350 B=EXP< <Sgy-Nvs*Sgf )/Nrun) 

3355 PRINT 

3360 PRINT USING "1 OX . ""Vapor-Shear Correlation: 

3365 PRINT USING "12X.""B = "".Z.3D":B 

3370 PRINT USING "12X.""n = . Z . 3D" : Nvs 

3375 END IF 

3380 ASSIGN GFilel TO * 

3385 ASSIGN «Fiie2 TO * 

3390 ASSIGN i«File3 TO * 

3395 SUBEND 
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3400 ! 

3405! THIS SUROUTINE CONVERTES THERMOCOUPLE VOLTAGE INTO TEMPERATURE 
3410! 

3415 SUB Tvsv(T) 

3420 COM /Ci/ C< 7 > 

3425 Siun'O . 

3430 FOR 1-0 TO 7 

3435 Sam* $u/n + C< I )*T I 

3440 NEXT I 

3445 T“Sum 

3450 SUBEND 

3455! 

3460! THIS FUNCTION CALCULATES PRANDTL NUMBER OF WATER IN THE 

3465! RANGE 15 TO 45 DEG C 

3470! 

3475 DEF FNPrw(T) 

3480 Y= 1 0 < 1 . 09976605- T*< 1 . 374932GE-2- T*< 3 . 968875E-5-3 . 45026E-7*T ) ) ) 

3485 RETURN Y 

3490 FNEND * 

3495! 

3500! THIS FUNCTION CALCULATES THERMAL CONDUCTIVITY OF WATER 

3505! IN THE RANGE OF 15 TO 105 DEG C 

3510! 

3515 DEF FNKu< T ) 

3520 Y* . 5625894+ T*< 2 .296 4546E-3-T«< 1 .509766E-5-4 . 0581652E-8*T ) ) 

3525 RETURN Y 
3530 FNEND 
3535* 

3540! THIS FUNCTION CALCULATES SPECIFIC HEAT OF WATER 

3545! IN THE RANGE 15 TO 45 DEG C 

3550! 

3555 DEF FNCpu.(T) 

3560 Y e ( 4 . 21 1 20858- T*(2 . 26826E- 3-T*< 4 . 4236 IE-5 + 2.71 428E -7*T >>>*1000 
3565 RETURN Y 
3570 FNEND 
3575! 

3580! THIS FUNCTION CALCULATES DENSITY OF WATER IN THE 

3585! RANGE 15 TO 105 DEG C 

3590! 

3595 DEF FNRhou>< T ) 

3600 Ro=999 . 52946+ T* < . 0 1 269- T* ( 5 . 48251 3E-3- T* 1 .234147E-5>> 

3605 RETURN Ro 
3610 FNEND 
3615! 

3620! THIS FUNCTION APPLIES CORRECTIONS TO THERMOCOUPLE READINGS 
3625! 

3630 DEF FNTemp(T.I) 

3635 DIM A(14).B(14> 

3640 DATA 0,640533,0.573054.0.593101 . 0 . 57298 . 0 . 56228 . 0 . 567384 , 0 . 569577 

3645 DATA 0,553951 ,0.552008.0.566955.0,520998.0.522661 .0.531008.0.560788,0.5524 

05 

3650 DATA 1 1 . 8744 ,8 . 631 63 . 9 . 394 1 2 .8 . 570246 . 8 . 299436 .8 . 36677 . 8 . 04507 . 7 .459766 

3655 DATA 7 . 498928. 7 .9408 . 5 .87072 .5 . 391 556 . 6 . 1 3399 . 6 . 48586 ,6 . 326224 

3660 READ A<*> .BH 

3665 IF I<15 THEN 

3670 T = T-(A( I ) - B < I>* . 00 1 *T ) 

3675 ELSE 
3680 T = T- .5 
3685 END IF 
3690 RETURN T 
3695 FNEND 
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3700! 

3705! THIS FUNCTION COMPUTES THE SPECIFIC VOLUME OF STEAM 
3710! 

3715 DEF FNVvst(Tt) 

3720 P-FNPvst < Tt ) 

3725 T-T t+273 . 1 5 

3730 X* 1 500/T 

3735 F 1 * 1 /< 1 +T* 1 .E-4> 

3740 F2«<1-EXP<-X>>'2.5*EXP<X)/X .5 

3745 B« .001 5*F1- .000942-F2- .0004882- X 
3750 K«2*P/<461 .52*T> 

3755 V* < 1 ♦< M2*B*K) * . 5 > /K 
3760 RETURN V 
3765 FNEND 
3770! 

3775! THIS FUNCTION CONVERTS THE VOLTAGE READING OF THE PRESSURE 

3780! TRANSDUCER INTO PRESSURE IN MM HG 

3785! 

3790 DEF FNPvsv(V) 

3795 Y*1 . 1103462+ 163. 3641 3*V 

3800 RETURN Y 
3805 FNEND 
3810! 

3815! THIS FUNCTION CALCULATES THE SATURATION TEMPERATURE OF STEAM AS A FUNCTION 

3820! OF PRESSURE 

3825! 

3830 DEF FNTvsp(P) 

3835 T u*110 
3840 T 1 =80 
3845 Ta=(Tu+Tl >* .5 
3850 Pc=FNPvst < Ta) 

3855 IF ABS((P-Pc)/P)>.0001 THEN 

3860 IF Pc<P THEN Tl = Ta 

3865 IF Pc>P THEN Tu=Ta 

3870 GOTO 3845 

3875 END IF 

3880 RETURN Ta 

3885 FNEND 

3890! 

3895! THIS FUNCTION COMPUTES THE VISCOSITY OF WATER 
3900! 

3905 DEF FNMuu(T) 

3910 A = 247 . 8/< T + 1 33 . 1 5 ) 

3915 Mu-2 . 4E-5* 1 0 'A 
3920 RETURN Mu 
3925 FNEND 
3930! 

3935! THIS FUNCTION COMPUTES THE LATENT HEAT OF VAPORIZATION 
3940! 

3945 DEF FNHf g< T > 

3950 Hf g=2497 . 7389-T «<?.2074 + T*< 1 . 7079E- 3-2 . 8593E-6-T ) ) 

3955 RETURN Hfg 
3960 FNEND 
3965 • 

3970! THIS FUNCTION COMPUTES THE SATURATION PRESSURE 
3975 1 

3980 DEF FNPvst ( Ts team > 

3985 DIM K ( 8 > 

3990 DATA -7 . 691234564 . -26 . 08023696 .-168.1 706546 . 64 . 23285504 . -118. 9646226 
3995 DATA 4 . 1 6711 732 . 20 . 9750676 . 1 E9 . 6 
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4UUU KtHO M*) 

4005 T - ( Ts team +27 3 . 1 5> /647 . 3 

4010 Sujn-0 

4015 FOR N-0 TO 4 

4020 Sum*Sam+K(N)*< 1-T)*(N+1 ) 

4025 NEXT N 

4030 Br-Sum/<T«<1+K<5>«<1-T>+K<6>«<1-T>*2>>-<1-T>/<K<7>-<1-T 

4035 Pr-EXP(Br) 

4040 P-22120000«Pr 

4045 RETURN. P 
4050 FNEND 
4055! 

4060! THIS FUNCTION CALCULATES THE STEAM MASS FLOW RATE 
4065! 

4070 DEF FNMdot(Dp) 

4075 Mdot-4.3183E-3+5.6621E-4«Dp 

4080 RETURN Mdot 
4085 FNEND 



>~2+K<8> ) 



* 
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APPENDIX D 



SEIDER-TATE CONSTANT CALCULATION 



The following solution procedure used in calculating 
the Sieder-Tate constant, by using a modified Wilson plot, 
is presented in this appendix: 



1. The overall heat- transfer coefficient is defined 
as : 1 1 D 



— = — + R w + — - — 
U h w D. h. 

o o 11 



(D.l) 



2. The Nusselt number is defined as: 





Nu 


= h. 

l 


D i / k c 




or 


Nu 


= c. 

l 


Re 0.8 pr l/3 (y c /y w ) 0 * 14 


(D.2) 


or 


Nu 


= c. 
1 


; where 


(D. 2a) 




Q. = 


Re° 


.8 D 1/3 , , .0.14 

Pr ^c /y w } 





By using the Nusselt equation: 

-,1/4 



h = 0.725 
o 



rk/ P f 2 h fg g 



p r D (T -T ) 
M f o v sat w 7 



and q = h (T -T ) 
M o v sat w 7 



(D.3) 

(D.4) 



and upon combining equations (D.3) and (D.4) 

1/3 



h = 0.651 

o 



3 2 , 

k f P f h fg 

P f D q q 



(D.5) 
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5. 



NOTE 



or h 



= 0.651 v 1 / 3 ; where 

i 3 2 , 

kf p f h fg g 



(D. 5a) 



v = 



Uf D 0 q 

Substituting equations (D.2a) and (D.5a) into 
equation (D.l), and rearranging gives: 

1 i/T 1 D v 1/3 1 

( R ) v i/3 = - + (D. 6) 

U w C. k n 0.0651 

o 1C 

or Z = mW+ 1/0.0651; where 



Z . . R ) v 1/3 , 

U o 

m = 1 / and 

W = D v 1/3 / k ft 
o c 

By plotting Z vs . W and computing the slope of the 
least-squares-fit line, and then taking the recip- 
rocal of this slope, the resulting value is the 
Sieder-Tate constant. 



T and T r must be determined iteratively, as 
w r 

described in Chapter IV, section E. 
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